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Background: Blood-brain equilibration rate constant (k e0 ) is derived from either pharmacokinetic and pharmacodynam- 
ic modeling (k e0 _ modd ) or a model-independent observed time to peak effect (k e0tpe> ^). Performance in bispectral index (BIS) 
prediction was compared between k e0 model and k e0 for microemulsion or long chain triglyceride (LCT) propofol. 
Methods: Time to peak effect (t peak , time to a maximally reduced BIS value) of microemulsion propofol after an intra- 
venous bolus (1 mg/kg) was measured in 100 patients (group A micro ). An observed t peak of 1.6 min for LCT propofol was 
obtained from an earlier study. Another 40 patients received a target controlled infusions of microemulsion propofol {k e0 _ modd 
= 0.187/min, group B micr0 = 20) or LCT propofol {k e0 moiA = 0.26/min, group B LCT = 20) and remifentanil. The fc^^'s in group 
B micr0 and B LCT were calculated using the observed tp eak value obtained from group A micro and 1.6 min, respectively. Effect- 
site concentrations of propofol were recalculated using the amounts of propofol infused over time and k e0 tpt ^s. Predicted 
BIS values calculated by sigmoid Emax equations with k e0 _ model and k e0 tpeak were compared with observed BIS values dur- 
ing induction and emergence for both formulations of propofol. 

Results: Observed t peak of microemulsion propofol was 1.68 min. The median performance errors of BIS in group B^,, 
were -1.83% (-24.8 to 18.9, k e0 _ model ) and -2.42% (-26.1 to 36.2, fc e0 tpe J, while 8.01% (-20.5 to 30.1, k e0 modd ) and 7.37% 
(-27.0 to 49.1, A: e0 tp eak ) in group B LCX . The median absolute performance errors of BIS in group B^^, were 11.87% (2.2- 
31.lW.dd) and 14.38% (-0.6 to 44.6, k e0 tpeak ), while 17.31% (5.54-36.0, k e0 _ model ) and 18.28% (-0.1 to 56.0, k e0 tpeak ) in 
group B LCT . 

Conclusions: The k e0 _ model showed better performance in BIS prediction than the k e0 tpeBk . (Korean J Anesthesiol 2013; 65: 
299-305) 
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Introduction 

Target-controlled infusion (TCI) systems were developed to 
administer intravenous agents during anesthesia. This is a logical 
approach for the fine titration of intravenous hypnotics and opi- 
oids and is based on the pharmacokinetic model that describes 
the disposition of the drug. For targeting effect-site concentra- 
tions, the TCI system should incorporate pharmacokinetic pa- 
rameters, such as the distribution volume and clearance, and the 
blood-brain equilibration rate constant {k e0 ). Pharmacokinetic 
parameters can be determined by assessing the time course of 
the drug. On the other hand, k e0 can be estimated using pharma- 
cokinetic/pharmacodynamic modeling of a single population 
(fc e0 model ) or by calculating model-independent time-to-peak ef- 
fects (fc e0 tpeak ) [1,2]. Although integrated pharmacokinetic/phar- 
macodynamic modeling of a single population is the standard 
method for obtaining k e0 , this is complicated by the need for 
blood samples and the frequent measurement of drug effects. As a 
result, the pharmacokinetic model for propofol that was reported 
by Marsh et al. [3], a popular model used in clinical situations, 
only estimates pharmacokinetic parameters without k e0 . 

The time of peak effect after bolus injection, t peak , could be 
used to link separate pharmacokinetic and pharmacodynamic 
studies [2]. The tp eak method is based on the assumption that the 
time required to reach the maximal effect-site concentration is 
identical to the time required to produce the maximal drug ef- 
fects. Minto et al. [2] reported that the method, which is an 
alternative approach that can be used to obtain k e0 , can also be 
used to obtain better approximations of the postulated "true" 
time course of propofol effects in comparison with the naive 
method. The observed t peak of propofol was 1.6 minutes in a 
previous study [1], which was then used to calculate k e0 using 
differential equations. The k e0 of long-chain triglyceride (LCT) 
propofol, which was calculated as 1.2152/min using the tp eak 
method, is included the Orchestra workstation (Fresenius Kabi 
AG, Germany) and was previously modified by Marsh et al. [4] . 

Microemulsion propofol (Aquafol; Daewon Pharmaceutical 
Co., Ltd, Seoul, Korea) was originally developed to overcome 
the drawbacks of LCT propofol [5]. Sequential pharmacoki- 
netic/pharmacodynamic modeling of microemulsion propofol 
was performed in healthy volunteers [6]. The comparison of 
appropriateness between /c e0 s derived from different methods in 
target effect-site concentration controlled infusion of propofol 
has yet to be fully elucidated. The aim of this study is to compare 
predictive performance in bispectral index (BIS, Aspect Medical 
Systems, Newton, MA, USA) prediction during target effect-site 
controlled infusion of propofol using different k e0 s in the same 
pharmacokinetic model. 



Materials and Methods 

Part 1: Patient population and study procedure used 
to determine the k e0 value of microemulsion propofol 
using the t peak method 

After obtaining approval from the Institutional Review Board 
of our institution and signed informed consent from all patients, 
a total of 100 patients diagnosed as American Society of Anes- 
thesiologist (ASA) physical status (PS) 1 or 2 who were sched- 
uled for stomach or colorectal surgery under general anesthesia 
were enrolled in this study. Patients with known or suspected 
cardiac, liver, or renal disease, significant obesity (> 120% of 
ideal weight), patients currently receiving sedative or analgesic 
drugs, and patients with known adverse effects to propofol were 
excluded. The patient characteristics are shown in Table 1. All 
patients were required to fast starting the midnight before the 
operation without receiving premedication. Once in the operat- 
ing room, patients were monitored using electrocardiography, 
pulse oximetry, end-tidal carbon dioxide partial pressure, 
noninvasive blood pressure monitoring (Datex-Ohmeda S/5; 
Planar Systems, Inc., Beaverton, USA), and BIS. Using RS232C 
cables, all data were continuously downloaded onto a personal 
computer until recovery from anaesthesia. An 18-gauge catheter 
was inserted into a large vein in the wrist. Forty milligrams of 
lidocaine were intravenously administered prior to propofol 
injection, then an intravenous bolus of 1 mg/kg microemulsion 
propofol was administered. After observing the maximal reduc- 
tion in BIS from baseline, propofol and remifentanil were ad- 
ministered by effect-site target-controlled infusion (Asan Pump, 
version 2.1.3; Bionet Co., Ltd., Seoul, Korea). 



Table 1. Patient Characteristics 

Part 2 

Part 1 

(n = 100) Microemulsion LCT propofol 
propofol (n = 20) (n = 20) 



Age (yr) 


52.0 (46.0-62.0) 


59.2 ± 10.2 


60.2 ± 10.2 


Height (cm) 


161.5 ±8.1 


159.7 ±9.6 


163.4 ±5.0 


Weight (kg) 


59.7 (54.6-67.6) 


59.2 ± 10.7 


63.8 ± 13.8 


Sex (M/F) 


52/48 


7/13 


15/5 


ASA PS (I/II) 


22/78 


3/17 


0/20 



Parts 1 and 2 of this study were designed to determine the k e0 value 
for microemulsion propofol using the t peak method. The predictive 
performance of BIS was determined during TCI using the 2 k e0 s 
values used in the same pharmacokinetic propofol model. Data are 
expressed as the mean ± standard deviation, median (range: 25-75%), 
or countswhere appropriate. LCT: long-chain triglyceride, ASA PS: 
American Society of Anesthesiologists physical status. 
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The differential equation describing effect-site concentration 
as a function of time t was presented as equation (1) 



dCe(t) 
dt 



K Q ■ Cp(t) - k e0 ■ Ce(i) = k e0 (Cp(t)- Ce(t)) 



(Eq.l) 



Equation 2 shows the Laplace transform of Equation (1). 

(Eq.2) 



sCe(s) = k e0 ■ Cp(s) - k e0 ■ Ce(s) 



Ce(s) = ^-Cp( S ) 

Equation 3 shows the inverse Laplace transform of the con- 
version of the time domain. 



Ce(t) = k e0 -e- k °°'*Cp(t) 



(Eq.3) 



Equation 4 shows the unit disposition function of the plasma 
concentration determined in a 3-compartment mammillary 
model. 



(Eq.4) 



Starting with Equation 3, the convolution of the input and 
disposition functions is shown in Equation 5. 

Ce(?) = k M jy^'-^Cpi^dT 

= C, •k e0 -e' K '-' jy-V^Vr = C t -k^-e 1 "'' jV'^'Vr 



--C.-k.-e 



--C,-k.„-e 



k„ n - X; 



e -e -e e ^ -e 



-I 



-I 



(Eq.5) 



As a result of Equation 5, the effect-site concentration of the 
3-compartment mammillary model can be summarized using 
Equation 6. 

Ce ( ? ) = X f"° '—. -e*« 4 ) (Eq.6), where Ci and Xi 



Table 2. Pharmacokinetic Parameters and k e0 s Values of the Noh and 
Modified Marsh Models 



Noh model 



Modified Marsh model 





21.0 


0.228 x WT 


V 2 (L) 


86.7 + (LBM / 47) 1L3 


0.473 x WT 


V,(L) 


757 


2.865 x WT 


CL (L/min) 


1.37 + 0.0091 x (WT - 64) 


0.027 x WT 


Q, (L/min) 


F: 0.855; M: 1.17 


0.026 x WT 


Ch (L/min) 


0.807 


0.010 x WT 


Ko (/min) 


0.187* 


1.2152 t 



F: female, M: male, LBM: lean body mass calculated using the James 
formula, WT: weight. *Estimated using sequential pharmacokinetic/ 
pharmacodynamic modeling of a single population. 1 Calculated based 
on the observed time-to-peak effects. 



are the coefficients and exponents of the pharmacokinetic mod- 
el, respectively. 

Because tp eak is the time at which effect-site concentration (Ce) 
reaches its maximal value, the derivative of Ce with respective to 
t should be 0 at t = tp eak . 

c<o=Zr^(V""*' w -Koe~* M *n=o (Eq.6) 

i=l k e0 

Thus, Equation 6 can be used to mathematically determine k e0 . 

Part 2: Patient population and study procedure used 
to compare BIS predictions during TCI based on the 
2 k e0 values determined in pharmacokinetic propofol 
models 

A total of 40 ASA PS 1/2 patients who were scheduled for 
stomach, breast, or colorectal surgery were randomly divided 
and assigned to receive microemulsion (n = 20) or LCT propo- 
fol (n = 20; Pofol, Jeil Pharmaceutical Co., Ltd, Seoul, Korea). 
The patient characteristics are shown in Table 1. All patients 
were required to fast starting the midnight before surgery with- 
out receiving any premedication. Once in the operating room, 
patients were monitored using electrocardiography, pulse ox- 
imetry, end-tidal carbon dioxide partial pressure, noninvasive 
blood pressure monitoring, and BIS. Using RS232C cables, these 
data were continuously downloaded onto a personal computer 
until recovery from anaesthesia. Each patient was preoxygenated 
with 100% oxygen via a facemask. Microemulsion propofol was 
administered via effect-site target-controlled infusion based on 
the Noh model, and a modified version of the Marsh model was 
used to administer LCT propofol to the other group of patients. 
The pharmacokinetic parameters and k e0 values of the Noh and 
modified Marsh models are shown in Table 2 [4,6]. Remifent- 
anil were administered via effect-site target-controlled infusion 
based on the Minto model [7]. Propofol and remifentanil were 
infused using TCI software, and the infusion histories of the 
patients were automatically recorded onto comma separated 
values (xsv) files. The pharmacokinetic parameters and k e0 val- 
ues shown in Table 2 were programmed into the Asan pump. 
Tracheal intubation was facilitated by administering 0.12 mg/ 
kg vecuronium. Patients were then ventilated with a mixture of 
oxygen and air (delivered at a 1 : 2 ratio), and the ventilation rate 
was adjusted to maintain the end-tidal carbon dioxide partial 
pressure between 35-45 mmHg. The target site-effective con- 
centration of propofol was adjusted to maintain BIS values < 60. 
Target site-effective concentrations of remifentanil were titrated 
to maintain stable hemodynamics (i.e., systolic blood pressure 
[SBP] > 80 mmHg and heart rate [HR] > 45 beats/min). If nec- 
essary, ephedrine or atropine was administered to maintain SBP 
> 80 mmHg and HR > 45 beats/min during anesthesia. Neuro- 
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muscular blockade was reversed by administering neostigmine 
and glycopyrrolate at the end of surgery. 

Part 3: Simulation based on the infusion propofol 
history described in Part 2 

The infusion history of propofol determined in Part 2 was 
used to simulate the effect-site concentration of propofol us- 
ing different k e0 values in same pharmacokinetic model. In part 
2, patients received microemulsion propofol using TCI with 
k e0 obtained by A: e0 estimation method. On the other hand, patients 
received LCT propofol using TCI with k e0 obtained by k^^y. 
method. For simulation, k e0 s of microemulsion and LCT pro- 
pofol was changed to different k e(l s obtained by k^j^ method 
for microemulsion propofol and A: e0 _ estimation method for LCT 
propofol. fc e0 estlation value for LCT propofol was obtained from 
the result of Billard et al. [8]. After each csv files was loaded in 
Asan Pump, k e0 s of propofol were changed to 0.26/min for LCT 
propofol and result of part 1 for microemulsion propofol [8], 
and simulations of individual patients were performed by Asan 
Pump. Model summary for TCI and simulation of microemul- 
sion and LCT propofol is shown in Table 3. 

Predictive performance 

In accordance with a previous method [9], the predictive per- 
formance of the TCI system was determined using the pharma- 
cokinetic parameters and the k e0 values for each type of propofol 
were characterized based on 2 parameters: bias and inaccuracy. 
BIS values (measured from the induction of propofol infusion to 
intubation) and recovery periods (measured from the discontin- 
uation of propofol infusion to extubation) were used to analyze 
the results of the present study. For each BIS value, the perfor- 
mance error (PE) of the ith patient was calculated as follows: 

measured;, - predicted,, 
J predicted^ n 

where predicted^ represents the predicted BIS value at the jth 
sampling point of the ith patient and measured^ is the observed 
BIS value. 



The sampling interval (tft^) for BIS measurement was 10 
seconds. The predicted value of BIS for each sample was calcu- 
lated using the inhibitory sigmoid E max model: 

Effect = E B +(E max -E 0 ) ^ ^ (Eq.8) 

where Effect is the BIS value, E 0 is the baseline BIS value when 
no drug is present, E max is the maximal possible drug effect on 
BIS, Ce is the predicted site-effective concentration of propofol 
determined in Part 3, Ce 50 is the site-effective concentration as- 
sociated with 50% of the maximal drug effect on BIS, and y is 
the steepness of the site-effective concentration vs the BIS rela- 
tionship. 

Various pharmacodynamic parameters for microemulsion 
and LCT propofol, such as E 0 , E max , Ce 5 „, and /, were determined 
in our previous studies [5,6]. The first performance measure, the 
bias (i.e., the direction and size of deviation from the target) of 
the ith individual, was calculated as the median performance er- 
ror (MDPEi) as follows: 



MDPE, = median PE, 



,J 



= 1, 



.N, 



(Eq.9) 



The inaccuracy of BIS prediction (i.e., the size of the typical 
miss) of the rth individual was calculated as the median absolute 
performance error (MDAPE;) as follows: 



MDAPE, = median \PE„ 



■ N, 



(Eq.lO) 



where AT f is the number of blood sampling points for the ith in- 
dividual. 

The pooled-data approach was used to obtain population 
estimates for bias and inaccuracy (fit4NM 3.5.1, Eun-Kyung Lee 
and Gyu-Jeong Noh, http://www.fit4nm.org/download, last ac- 
cessed: Oct 17, 2011). 

Statistical analysis 

All statistical analyses were conducted using R software (ver- 
sion 2.13.1; R Foundation for Statistical Computing, Vienna, 
Austria) or SigmaStat 3.5 for Windows (Systat Software, Inc., 
Chicago, IL). Data are expressed as the mean ± standard devia- 
tion (SD) for normally distributed continuous variables, the 



Table 3. Model Summary of the Target-controlled Infusion and Simulated Infusion of Microemulsion Propofol and Long-chain Triglyceride Propofol 

Microemulsion propofol (n = 20) LCT propofol (n = 20) 

Pharmacokinetic model of propofol 
Pharmacokinetic model of remifentanil 
k e0 for target-controlled infusion (min) 
Target effect-siteconcentration of propofol (ug/ml) 
k e0 for simulation based on infusion history 

LCT: long-chain triglyceride, A: e0 estiIlation : estimated from sequential pharmacokinetic/pharmacodynamic modeling of a single population, k t 
calculated using the observed time-to-peak effects. 



Noh model 
Minto model 



eO_estimation 

1.0-2.6 



0.187 



Time-to-peak effect: 1.68 min 



Modified Marsh model 
Minto model 
Kojpctt- 1-2152 
2.0-3.0 
Kn ahJ 0.26 /min 
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median (range: 25-75%) for nonnormally distributed continu- 
ous variables, and as counts and percentages for categorical 
variables. The relationship between age, time-to-peak effects, 
and BIS was analyzed using linear regression. All P values < 0.05 
were considered statistically significant. 

Results 

All patients enrolled in the study of obtaining k e0 of micro- 
emulsion propofol calculated by t peak method were included 
in order to analyze the time to peak effect. Time course of BIS 
values after an intravenous bolus injection of microemulsion 
propofol 1 mg/kg are shown in Fig. 1. Median (25-75%) value 
of the time to maximally reduced BIS from baseline value was 
1.68 (1.24-2.45) min. Also, relationship between age and time to 
peak effect in BIS are presented in Fig. 2, where, with increasing 
age, time to peak effect (tpj increased. 

A total of 1,027 BIS values (average value: 63.3 ± 16.9, sam- 
pling duration: 17.1 ± 3.5 min) from 20 patients receiving mi- 
croemulsion propofol and 1,878 BIS values (average value: 59.3 
± 17.7, sampling duration: 24.0 ± 7.2 min) from 20 patients re- 
ceiving LCT propofol were used to determine predictive perfor- 
mance in BIS prediction. Performance error (PE) and absolute 
performance error (APE) in BIS values for microemulsion pro- 
pofol and LCT propofol according to the different k e0 s were de- 
picted in Fig. 3. PE of ^estimation was significantly lower than that 
of fc^jp^ in microemulsion propofol (P = 0.013, Mann- Whitney 
rank sum test). APE of /c e(L( , stimation showed significant differences 
in both formulations, compared with k e0J[peilk (P < 0.001, Mann- 
Whitney rank sum test). Bias (median performance error) and 
inaccuracy (absolute performance error) in BIS prediction val- 
ues for microemulsion and LCT propofol in same pharmacoki- 



netic model with different k e0 s are demonstrated in Table 4. BIS 
values predicted by ^estimation method showed higher accuracy 
than those predicted by fc e0 tpeak method. 

Discussion 

The median t^ of microemulsion propofol was 1.68 min. 
The k e0 derived from sequential pharmacokinetic and pharma- 
codynamic modeling showed higher accuracy in BIS prediction 
than that of k e0 calculated by the time to peak effect in both pro- 
pofol formulations. 

The tp eak method should be used only for linking separate 
pharmacokinetic and pharmacodynamic studies for the same 
drug, as employed in an earlier study [1], where tp eak increased 
with age and this finding was consistent with our results. In 
other study [10], t peak was not affected by age, which may be ex- 
plained by the fact that distribution of age was relatively narrow 
(49 ± 16 yr). Age of patients enrolled in present study for obtain- 
ing k e0 of microemulsion propofol calculated by t^ak method was 
24-86 yr, which was similar to the previous study (25-81 yr) [1]. 
Overall, elderly patients showed prolonged onset of action on 
administration of drug [11,12]. Age was a significant covariate 
for k e0 in population pharmacodynamic modeling of remifent- 
anil, where t 2n k e0 increase with increasing age [7,13]. These facts 
may be explain the relationship between age and time to peak 
effect, even though effect of age on k eQ of propofol using popula- 
tion pharmacodynamic analysis has not yet been elucidated. 

In Noh model, lean body mass, weight and sex were signifi- 
cant covariates for V 2 , CL and Q p respectively, resulting in differ- 
ent hybrid constants according to the individual patients. Hence, 
k e0 s of microemulsion propofol calculated by t peak method have 
different values when height, weight and sex are not identical 




012345678 
Time (min) 

Fig. 1. Raw recordings of the bispectral index (BIS) in 100 patients who 
received an intravenous injection of 1 mg/kg microemulsion propofol. 
The thick lines and black dashed lines indicate the mean and observed 
BIS values of individual patients, respectively. 
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Fig. 2. Relationship between age and time-to-peak effects on the 
bispectral index (BIS). The solid lines and black dashed lines indicate the 
regression line and 95% confidence intervals, respectively. Time-to-peak 
effect = 0.03 x age (years) + 0.74 (simple linear regression; P = 0.03). 
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Fig. 3. Performance error and absolute performance error of the bispectral index (BIS) values for microemulsion propofol and long-chain (LCT) 
triglyceride propofol according to 2 k e0 s values. k e0 _ madA : estimated from sequential pharmacokinetic/pharmacodynamic modeling of a single 
population, k ea _^ eak : calculated using the observed time-to-peak effects. The values of the k e0 _ moid and k e0 _ tpal models were 0.187/min and 1.68 min 
(time-to-peak effects) for microemulsion propofol and 0.26/min and 1.2152/min for LCT propofol, respectively. *P < 0.05 vs k e0J?c2k . Performance 
error (PE) was calculated as follows: 

measured BIS tj - predicted BIS tj 
' predicted BISjj 

Predicted BIS was calculated using the following equations. For microemulsion propofol, 

effect-site concentration 3 05 



Predicted BIS = 87.8 + (22.2 - 87.8)- 



For LCT propofol, 
Predicted BIS = 83.0 + (27.5 



■ 83.0)- 



1.58" + effect-site concentration 



effect-site concentration 4 3 



1.95 + effect-site concentration 

The pharmacodynamic parameters for microemulsion and LCT propofol that were used to calculate the predicted BIS values were obtained in 
previous studies. Empty circles (O) and red horizontal lines indicate the performance error of each BIS measurement and the median values, 
respectively. 



Table 4. Comparison of the Biases and Inaccuracies of Bispectral Index Prediction 





Microemulsion propofol (n = 20) 


LCT propofol (n = 20) 




k k 

fv e0_model fv e0_tpeak 


^e0_model ^e0_tpeak 


Bias (MDPE), % (95% CI) 
Inaccuracy (MDAPE), % (95% CI) 


-1.83 (-24.81 to 18.95) -2.42 (-26.10 to 36.18) 
11.87 (2.16 to 31.13) 14.38 (-0.61 to 44.60) 


8.01 (-20.55 to 30.06) 7.37 (-27.03 to 49.10) 
17.31 (5.54 to 35.98) 18.28 (-0.10 to 56.04) 



LCT: long chain triglyceride, MDPE: median performance error, MDAPE: median absolute performance error, k e0 mode i~ estimation from sequential 
pharmacokinetic-pharmacodynamic modeling in a single population, k e0Jpal : calculation using the observed time to peak effect. The value of k e0 _ moiA 
and k M X^ for microemulsion propofol and LCT propofol were 0.19 /min and 1.68 min (time to peak effect) for microemulsion propofol and 0.26 /min 
and 1.21 /min for LCT propofol, respectively. 
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each other. On the other hand, in modified Marsh model, weight 
was a significant covariate for all pharmacokinetic parameters in 
volume and clearance domain, resulting in same hybrid constant 
due to reducing a fraction with weight. The k e0 of LCT propofol 
calculated by tp eak method is identical (1.2152 /min) in modified 
Marsh model, irrespective of patients' demographics. 

In general, predictive performance of a TCI system was clini- 
cally acceptable when inaccuracy is less than 20-30 % [9]. Phar- 
macodynamic predictions as described by BIS exhibited smaller 
bias and inaccuracy in present study. Theoretically, because 
effect-site concentration per se should be statically related to ef- 
fect, target effect-site controlled infusion using Noh or modified 
Marsh models can be available to predict BIS approximately, ir- 
respective of method for obtaining k e0 . 

Biases and inaccuracies of a TCI system may be small- 
est when it incorporates pharmacokinetic parameters and k e0 
obtained from a single population. Potential causes of perfor- 
mance errors in the predictive performance of a TCI system 
may include model misspecifications and mechanical errors 
of the infusion pump, such as trumpet curves [14]. A trumpet 



curve could be more likely to develop on long time of infusion. 
In this study, effect of trumpet curve can be excluded from in- 
terpretation of result, because predictive performance of BIS 
prediction was analyzed in induction and recovery periods. 
Model misspecification of pharmacodynamic model shows a 
tendency of being smaller in the sequential pharmacokinetic 
and pharmacodynamic modeling, compared with simultane- 
ous pharmacokinetic and pharmacodynamic modeling. Model 
misspecification of pharmacokinetic model can convert to that 
of pharmacodynamic model in simultaneous pharmacokinetic 
and pharmacodynamic modeling. Pharmacokinetics of micro- 
emulsion propofol were well described by mammillary three 
compartment model and the results of typical pharmacokinetic 
parameter were used to characterize the relationship concentra- 
tion and effect using inhibitory sigmoid Emax model. 

In summary, the median time to maximally decreased BIS 
from baseline value after an intravenous bolus of microemulsion 
propofol 1 mg/kg was 1.68 min. BIS values predicted by fc c0 _ estimation 
method showed higher accuracy than those predicted by k e0 _ tpeilk 
method. 
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